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Abstract 

CO ' For each positive number a we study the analog v a in free probability of the 

classical Gamma distribution with parameter a. We prove that v a is absolutely 
\ continuous and establish the main properties of the density, including analyticity 

' and unimodality. We study further the asymptotic behavior of v a as a I 0. 

<D 

m ■ 1 Introduction 



In this paper we study the free Gamma distributions, i.e., the images of the classical 
Gamma distributions under the bijection between the classes of infinitely divisible mea- 
sures in classical and free probability, respectively, introduced by Bercovici and Pata (cf. 



|BP99] and [B NT02a] ). More precisely, for any positive number a the free gamma distri- 
bution v a with parameter a is defined as A(ji a ), where A is the Bercovici-Pata bijection 
(see Section F2]) and fi a is the classical Gamma distribution with parameter a, i.e., 

^ a {B) = -^-f r-V'dt (1.1) 

m . 

r**~ . for any Borel set B in 

The classical Gamma distributions form perhaps the simplest class of selfdecomposable 
measures on R which are not stable (see Section [2]). Since A preserves the notions of 
stability and selfdecomposability (see [BP 92] and [BNT02aj), the measures u a are thus of 
interest as (the simplest?) examples of non-stable selfdecomposable measures with respect 
to free (additive) convolution. Of particular interest is the free ^-distribution A(xf), 
^ ! which (up to scaling by 2) equals the measure v\ji- Apart from the general importance of 
the ^-distribution in classical probability, this is mainly due to the fact that the square 
of the semi-circle distribution (the analog of the Gaussian distribution in free probability) 
equals the free Poisson distribution (the image of the classical Poisson law by A) as 
observed e.g. in [VDN92J. Since A is injective, the relationship between the Gaussian and 
the ^-distribution thus breaks down in free probability, and from that point of view it is 
of some interest to identify further the measure A(x^). 

In an appendix to the paper [BP99J P. Biane studied the freely stable distributions and 
established their absolute continuity (with respect to Lebesgue measure) as well as the 
main features of their densities; in particular analyticity and unimodality. Applying the 
same method as Biane (based on Stieltjes inversion) we establish in the present paper that 
the free Gamma distributions v a are absolutely continuous with analytic densities, and 
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that they have supports in the form [s a , oo) for some strictly positive number s a , which 
increases (strictly) with a and tends to and oo as a goes to and oo, respectively. We 
derive an (implicit) expression for the density f a of v a in the form: 

fa(Pa(x)) = - 2 l aiX ] V2 (X G [-C*, OO)), (1.2) 

where P a is a strictly increasing function given by 

f°° t 2 e"' 

P Q (x)=2x + a-a —5— -—r dt, (x G [—c a , 00)), 

J {x-ty + v a {xy 

and c a is a positive constant such that P a {— c a ) = s a . Moreover v a : R — > R is a function 
essentially defined by the condition: 



x 



+ \v a {x)) [l + aG IMl (x + iv a (x))] G R, 



where denotes the Cauchy (or Stieltjes) transform of a probability measure /1011K 
(see formula (I2.3P below). This condition emerges naturally from the method of Stieltjes 
inversion in combination with the key formula: 

G Va {zO- + aG„ 1 {z))) = -, (1.3) 



which holds for all z in C + satisfying that z(l + aG fll (z)) G C + . The passage from (jl.3p 
to fll.2p via Stieltjes inversion depends heavily on a fundamental result of Bercovici and 
Voiculescu, which we state in Lemma [2.11 below for the readers convenience. 

By careful studies of the functions v a , P a and the right hand side of (II. 2p . we derive 
some main features of the density f a , e.g. analyticity, unimodality and the asymptotic 
behavior: 

l im J^L = ae ^ and lim ( f a ^) \ = Vg . (L4) 

In particular it follows that u a has moments of all orders, which is in concordance with 
the results of Benaych- George in (BG06J. 

We study also the asymptotic behavior of v a as a j. 0, and we prove that the measures 
—v a converge to the measure x~ l e~ x l^_ OQ )(x) dx in moments and in the sense of point- wise 
convergence of the densities: 

lim Ml = £-i e -« (£G(0,oc)). 

aiO a 

The remainder of the paper is organized as follows: In Section [2] we collect background 
material on infinite divisibility, The Bercovici-Pata bijection and Stieltjes inversion. In 
Section [3] we establish absolute continuity of v a and prove the expression (II .2p for the 
density. In Section H] we establish the asymptotic behavior ( II .4p and study how the 
quantities c a and s a vary as functions of a. In Section [5] we prove that v a is unimodal, 
and in the final Section [6] we study the asymptotic behavior of v a as a \. 0. The main 
results in Sections [3J1S] depend in part on some basic properties of the functions v a and P a , 
the proofs of which are (not surprisingly) rather technical. In order to maintain a steady 
flow in the paper, these proofs are deferred to Appendix |A] at the end of the paper. 
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2 Background 

2.1 Classical and free infinite divisibility 

A (Borel-) probability measure p on R is called infinitely divisible, if there exists, for each 
positive integer n, a probability measure p n on R, such that 

P = * [l n * ■ ■ ■ * jjy (2.1) 

n terms 

where * denotes the usual convolution of probability measures (based on classical inde- 
pendence). We denote by JT>(*) the class of all such measures on R. 

We recall that a probability measure p on R is infinitely divisible, if and only if its 
characteristic function (or Fourier transform) p has the Levy-Khintchine representation: 

fi{u) = exp irju - \au 2 + / (e iut - 1 - iutl { _ hl] (t)) p(dt)l , (u G R), (2.2) 

where r] is a real constant, a is a non-negative constant and p is a Levy measure on R, 
meaning that 

p({0}) = 0, and / min{l,t 2 } p(dt) < oo. 
Jr 

The parameters a, p and 77 are uniquely determined by p and the triplet (a, p, 77) is called 
the characteristic triplet for p. 

For two probability measures p and z/ on R, the free convolution // ffl v is defined as 
the distribution of x + y, where x and ?/ are freely independent (possibly unbounded) self- 
adjoint operators on a Hilbert space with spectral distribution p and z/, respectively (see 
|BV93j for further details). The class JD(EEI) of infinitely divisible probability measures 
with respect to free convolution ffl is defined by replacing classical convolution * by free 
convolution ffl in (12. ip . 

For a (Borel-) probability measure p on R with support supp(/i), the Cauchy (or 
Stieltjes) transform is the mapping G M : C \ supp(/i) — > C defined by: 

G,{z) = [ — + Mdt), (zeC\ supp(/i)). (2.3) 

The /ree cumulant transform C M of /1 is then given by 

e,(z) = zG^\z)-l (2.4) 

for all z in a certain region i? of C~ (the lower half complex plane), where the (right) 
inverse G^ 1 ^ of G^ is well-defined. Specifically R may be chosen in the form: 

R = {z e C~ I ~ G A^m}, where A^at = G C + | |Re(«)| < r]\m(z), \m(z) > M} 

for suitable positive numbers 77 and M. It was proved in |BV93j (see also }Ma92] and 
|Vo86] ) that C M constitutes the free analog of log/ M in the sense that it linearizes free 
convolution: 

G^mu(z) = e M (*) + Q v (z) 

for all probability measures p and v on R and all z in a region where all three transforms 
are defined. The results in |BV93j are presented in terms of a variant, (p^, of C M , which 
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is often referred to as the Voiculescu transform, and which is again a variant of the R- 
transform introduced in |Vo86j. The relationship is the following: 



= Mi) = «e M (i) (2.5) 

for all z in a region as above. In [B V93j it was proved additionally that p G JB(EB), 
if and only if there exists a in [0, oo), 77 in R and a Levy measure p, such that C M has the 
free Levy-Khintchine representation: 

Q^z) =r] z + az 2 + [ (—^— - 1 - tel t -i,i](*)) P(dt). (2.6) 

(cf. also [ BNT02b| ). In particular it follows for p in JD(EB) that C M can be extended to 
an analytic map (also denoted C M ) defined on all of C~. The triplet (a,p,T]) is uniquely 
determined and is called the free characteristic triplet for p. 

It was proved in |BV93| Proposition 5.12] that any measure v in JD(ffl) has at most 
one atom. In fact the proof of that proposition reveals that an atom a for v is necessarily 
equal to the non-tangential limit of (p v {z) as z — > 0, z G C + . We say that a function 
u: C + — > C has a non-tangential limit £ at 0, if for any positive number 5 we have that 

£ = lim u(z), where A s = {z G C + | lm(z) > 5\Re(z)\}. (2.7) 

In order to derive non-tangential limits, the following lemma (Lemma 5.11 in |BV93j ) is 
extremely useful: 



2.1 Lemma ([BV93J). Let u: C + — > C + be an analytic function, and let T be a curve 



in C + which approaches non-tangentially. 

If lim z ^o,zer u (z) = £, then lim^^o.^eAa u(z) = £ for any positive number 5, i.e., u has 
non-tangential limit £ at 0. 



2.2 The Bercovici-Pata bijection 

In |BP99] Bercovici and Pata introduced a bijection A between the two classes JD(*) and 
JD(EB), which may formally be defined as the mapping sending a measure \x from JD(*) 
with characteristic triplet (a, p, 77) onto the measure A(/i) in JD(EB) with free characteristic 
triplet (a, p, r\). It is then obvious that A is a bijection, and it turns out that A further 
enjoys the following properties (see [BP99J and |BNT02a] ): 

(a) If fi U fjL 2 G JD(*), then A(p x * p 2 ) = A(jii) ffl A(/i 2 ). 

(b) If p G JT>(*) and c G R, then A(D c p) = D c A(p), where e.g. -D c p is the transformation 
of p by the mapping x 1— > cx : R — > R. 

(c) For any constant c in R we have A(S C ) = 5 C , where 5 C denotes Dirac measure at c. 

(d) A is a homeomorphism with respect to weak convergence. 

Most of these properties can be established rather easily from the following convenient 
formula: 

POO 

Q AM (iz)= \ogfi(zx)e- x dx, (z G (—00, 0), ji G J2)(*)), (2.8) 
Jo 
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which was derived in |BNT04] . The properties (a)-(c) imply that A preserves e.g. the 
classes of stable and selfdecomposable measures. Specifically, let CP denote the class of 
all (Borel-) probability measures on M, and recall then that a measure \x from 7 is called 
stable, if it satisfies the condition: 

Va, a' > 3a" > 3/3 G R: D a p * D a ,p = D^p * 5 P . (2.9) 

Recall also that p is selfdecomposable, if 

Vc G (0, 1) 3p c E 7: p = D c p * p c . (2.10) 

Denoting by S(*) and £>(*) the classes of stable and selfdecomposable measures, respec- 
tively, it is well-known (see e.g. |Sat99] ) that §(*) C £(*) C JD(*). The classes S(ffl) 
and £(ffl) are defined be replacing classical convolution * by free convolution EB in f)2.9jl - 
([2TT0j) above. It was shown in jBV93j and jBNT02aj that S(ffl) C £(ffl) C JD(ffl). By 
application of properties (a)-(c) above, it follows then easily that 

A(S(*)) = S(ffl), and A(£(*)) = £(ffl) (2.11) 

(see |BV93j and [BNT02a] ). The measures in §(*) may alternatively by characterized as 
those measures in JD(*) whose Levy measure has the form 

P (dt) = (c_|t|- 1 - a -l ( _ oo , )(t) +C + r 1 - a + 1(0,00) (t)) dt 

for suitable numbers c + , c_ in [0, oo) and a + , a_ in (0, 2). Similarly £(*) may be character- 
ized as the class of measures in JD(*) with Levy measures in the form: p(dt) = \t\~ l k(t) dt, 
where k: M\{0} — » K is increasing on (— oo, 0) and decreasing on (0, oo). By the definition 
of A and (12. lip we have the exact same characterizations of the measures in S(ffl) and 
£(EB), respectively, if we let the term "Levy measure" refer to the free Levy-Khinthcine 
representation (12.61) rather than the classical one (12. 2|) . 

For any positive number a, the classical Gamma distribution fi a with parameter a (cf. 
(II. ip ) has Levy measure 

p a (dt) = Qft _1 e~*l(o )00 )(t) dt, 

and thus \i a G \ §(*). The corresponding free Gamma distribution, v a = A(// a ), 
satisfies accordingly that u a G £(ffl) \ S(ffl). As mentioned in the introduction, the 
purpose of the present paper is to disclose the main features of u a for any a in (0, oo). 

2.3 Stieltjes inversion. 

Let fi be a (Borel-) probability measure on R, and consider its cumulative distribution 
function: 

F^t) =M((-oo,t]), (t GR), 
as well as its Lebesgue decomposition: 

/j = p + cr, 

where the measures p and a are, respectively, absolutely continuous and singular with 
respect to Lebesgue measure A on K. It follows from De la Valle Poussin's Theorem (see 
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[Sak37j. Theorem IV. 9. 6]) that p and a may be identified with the restrictions of \x to the 
sets 

D 1 = {x E R I \mi h ^ F ^ x+h) ~ F ^ x) exists in R} 

and 

D O0 = {xER | lim^ ^±^M = oo}, 
respectively. In addition we have that (see e.g. Theorem 3.23 and Proposition 3.31 in 

HSU) 

A(R \ D x ) = 0, and p(dt) = i^(t)l Dl (t) dt, 

where, for any t in D\, F' (t) denotes the derivative of F^ at t. 

Consider now additionally the Cauchy (or Stieltjes) transform defined in (12. 3p . It 
follows then from general theory of Poisson-Stieltjes integrals (see [Do62j) that 

F 'Ax) = lim \m(GJx + iy)) for all x in Dx, 

71 yiO 



and that 



lim I lm(G M (x + iy)) I = oo for all x in D^. 



In particular we may conclude that the singular part a of p is concentrated on the set 

[x e R | lim^o \Gp(x +iy)\ = oo} 
(see also Chapter XIII in [RSTS] ). 



3 Absolute continuity of u a 



In this section we establish absolute continuity of the free Gamma distributions u a , a > 0, 
and prove the formula (11.21) for the densities. Our starting point is the derivation of the 
formula (11.31) . and we introduce for that purpose the function H a : C \ [0, oo) —> C given 
by 

f°° e~* f°° te~* 

H a {z) = z + zaG w (z) = z + za / dt = z + a + a / dt (3.1) 

Jo z — t J z — t 

for z in C \ [0, oo). By differentiation under the integral sign, note that H a is analytic on 
C \ [0, oo) with derivatives given by 



H' a (z) = l-a — dt, («6C\[0,oo)), 

/°° fe _ * 
, _ t)k+1 dt, (z6C\[0,oo), fee {2,3,4,...}). 



In the following we consider in addition the function F: C \ [0, oo) — > (0, oo) given by 

F(x + iy)= -df= / 7 ^r—jdt (3.3) 



\x + iy — t\ 2 J (x — t) 2 + y 2 
for all x,y G M such that x + iy E C \ [0, oo). 
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3.1 Lemma. Let a be a positive number. 

(i) There exists a unique positive real number c a such that 

- = F(-c a ) = / dt. 
a Jo [pa + t) 2 

The number c a increases with a, and satisfies that 

lim c a = 0, and lim c a = oo. 

(ii) There is a function v a : R — > [0, oo), such that 

{z e C + | Hjz) e C + } = {x + iy | x, y e R, y > w Q (x)} 
Tie function v a is given by 

v a (x) = 0, if x G (— oo, — c a ], 

F(x + ifafx)) = — , if x E (—c a , oo). 
a 

(iii) For all x in R we have that H a (x + ifo,(x)) G R. 

(iv) Tie function v a satishes that v a (x) > for all x in (— c a , oo). 




The graphs of the functions K1/2, ^x, «2 and vio- 

Proof of Lemma \3.1i (i) The function x t— >• F(—x) = J °° ,*^ 2 dt is clearly strictly 
decreasing and (by dominated convergence) continuous on (0, oo). Moreover, by monotone 
convergence, 

limF(— x) = oo, and lim F(—x) = 0. 

Hence, there exists a unique number c a in (0, oo) such that F(— c a ) = -. The last 
assertions in (i) are immediate from the last equality in ( 13. 4p . 
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(3.4) 

(3.5) 

(3.6) 
(3.7) 



(ii) For x, y in M such that x + iy G C \ [0, oo) we find from formula (13. ip that 

/ r°° f e -t \ 

lm(fl' a (ar + i2/)) = 2/ + alm( / — — -dt 

V J x + iy-t ) 

= y-ay / — - — - dt = y(l - aF(i + iy)). 

Jo i x - 1) + V 



For fixed i in 1 the function y i— )■ F(x + iy) is clearly strictly decreasing on (0, oo), and 
F(x + iy) — > as y — > oo. Moreover, by monotone convergence, 

i t-,/ \ I °°) if x > 

limFfa; + iy) = < 

2/40 v yi \F(x), ifx<0. 

Thus, if x > —c a , then lim y ^ F(x + iy) > F(—c a ) = -, and there exists a unique y x in 
(0, oo) such that F(x + iy x ) = -. Thus, if we put v a (x) = y x , then aF(x + iy) < 1 and 
hence lm(iJ Q (x + iy)) > for all y in (v a (x), oo). Similarly lm(if a (x + iy)) < for y in 
(0,f a (x)). 

If a; < — c a , then for all ?/ in (0, oo) we have that F(x + iy) < F(x) < F(—c a ) = -, 
and hence that \m(H a (x + iy)) > 0. Thus, if we put v a (x) = for x in (— oo, — c a ], it 
follows altogether that v a satisfies (13.51) . and that v a is given by (I3.6l) - (l3.7p . 

(iii) follows immediately from (13. 8p in combination with (I3.6p - (l3.7p . and (iv) is a 
consequence of the way v a was defined in the proof of (ii). ■ 

In continuation of Lemma [3.11 we introduce next the the following notation: 

S a = {x + iv a (x) \ x e R} (3.9) 

S' Q = {X + iV a {x) I X E [-Ca, Oo)} = Sa \ (-OO, ~C a ). (3.10) 

Sa — { x + iy I x i y G K, and y > f a (x)}. (3-H) 

Note in particular that ^ S a (since v a (0) > 0), and that S a , C C \ [0, oo). 

3.2 Proposition. Let a be a positive number, and consider the free Gamma distribu- 
tion v a with parameter a. Consider further the classical Gamma distribution /ii with 
parameter 1 (cf. (11. ip ). We then have (cf. formulae ( 12. 3p and ( 12. 4 p j 

(i) e i , Q (i) = «G' Ml (^)foraII^inC + . 

(ii) G Ua (H a (z)) = lforallzinSt. 

Proof, (i) The classical Gamma distribution fi a has characteristic function 

-t 



fJ, a [U) 



/ r°° e^' \ 

exp (^a J (e mt - l) — dtj , (« 6 1), (3.12) 

(see e.g. [Sat99j Example 8.10]). By formula (I2.8P and Fubinis Theorem it follows then 
for any u in (— oo, 0) that 

Q Ua {iu)= log jl a (ux)e- x dx = a (/ (e iuxt - I) — dt)e~ x dx 
Jo Jo ^ Jo t J 

f e"V 1 \ . Z" 00 e" 4 , 
a / — ■ 1 \ dt = am / ■ — dt. 



t V 1 — iut J ./n 1 — iw£ 
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Setting u = — ^, we find for any y in (0, oo) that 

roc Q -t r-co Q -t 

G u (^-) = a — ; t- dt = a dt = aG u , (i-u), 

and by analytic continuation we conclude that G Ua (l/z) = aG^z) for all z in C + . 
(ii) Recall from the definition of Q Ua (see Subsection 12. 1 p that 

e„ (i) = ^G^a) - 1 

for all 2 in a region of the form A v ^m = {z G C + | [ Re(^r) | < r]\m(z), \m(z) > M} for 
suitable positive numbers n and M. Taking (i) into account, we find that 

Gi a 1) (l) = z + zaG^ 1 (z) = H a (z), and hence ± = G Va (H a (z)) (3.13) 

for all z in A^m- Since G Va and H a are analytic on C + , it follows from Lemma f3. llf ii) 
and analytic continuation that the latter equation in (I3.13P holds for all z in 9J- This 
completes the proof. ■ 

In order to combine Proposition I3.2( ii) with the method of Stieltjes inversion (and 
Lemma I2.ip . we need some further preparations, which are presented in the series of 
lemmas to follow. 

3.3 Lemma. For any positive number a we have that 

H' a {-c a ) = 0, and H' a {z) ^ for all z in $ a \ {-c a }. 

In fact, 

d 

Re(H'(x + ivJx))) = -avJx)—F(x + ivJx)) > 

oy 

for all x in (—c a , oo). 

Proof. Note first that by (|52]l - (pPJ) we have that 

(c7+ ty 

For z = x + iy in C\ [0, oo) we find next, by application of the Cauchy-Riemann equations 
and (JE2D, that 

Re(i£(*)) = JW«W) = ^™( H *(z)) 

= —(y(l - aF(z)) = (1 - aF{z)) - ay-^F(z). 

For any x in (—c a , oo) it follows thus from (13. 7p that 

9 

Re(H' a (x + iv a (x)) = - av a (x)—F(x + iv a (x)). 

The proof if concluded by noting that differentiation with respect to y in (I3.3P leads to 

d w . f°° te-* 

— F(x + ly) = -2y / - -r— — — dt, 

dy J {{x-ty + y 2 ) 2 



H' a (-c a ) = 1 - a I -_ ; dt = 1 - aF(-c a ) = 0. 



o 



where the right hand side is strictly negative whenever y > 0. 
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In the following lemma we collect some further properties of the function v a , that will 
be needed in various parts of the remainder of the paper. We defer the rather technical 
proof to Appendix |A] 

3.4 Lemma. Let a be a positive number and consider the function v a : R — > [0, oo) given 
by fl3.(jj) -( ]5TT|) . Then v a has the following properties: 

(i) v a is continuous on R and analytic od! \ {— c a }. 

/• • \ r v a (x) 

[11) lim — — = air. 

x— >oo XG x 

(hi) For any positive numbers 5, 7 there exists a positive number a such that 

v a (x) 



sup 

x-S[<5,oo) 



— irxe x 



a 



< 7, whenever a G (0, cto]. 



3.5 Lemma. Consider a fixed positive number a. 

(i) For any z in (—00, — c a ) we have that z + it G St f° r a ^ positive t, and that 
H a (z + it) — > H a (z) G R non-tangentially (from C + to R), as t \. 0. 

(ii) For any point z in S' a \ {— c a } there exists a vector ^ z in C and a number e z in 
(0, 00), such that 

(a) z + t^ z G Sq for all t in (0, e z ). 

(b) H a (z + t^ z ) -)■ H a (z) G R non-tangentially (from C + to R), as 1 1 0. 

Proof, (i) Assume that 2 G (— 00,— c Q ). According to Lemma 13.1 ( ii) we have that 
z + it G So and hence + it) G C + for all positive t. It remains then to show that 

\m(±H a (z + it)) ^ at t = 0. Using (Q we find that 



^H a (z + it) = iH' a {z + it)=i(l-a J 



se 

ds 



(z + it- s) 2 

(z — s — it) 2 se~ 
{{z-s) 2 + t 2 ) 2 



ds 



and hence at t = we have that 



/ d \ f°° se~ s 

\m( M H a (z + it))=l-aJ o ¥ ^y 2 ds>0, 

since z < —c a (cf. f)3.4p ). 

(ii) Assume that z = x + iv a (x) for some x in (— c a , 00). We choose then a real number 
r, such that 

(1) r > v' a (x) (cf. Lemma l3T4T i)). 

(2) The vector (l,r) is not perpendicular to the vector (I m(H' a (z)), Re(H' a (z))) in M. 2 (cf. 
Lemma 13.31) . 
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We then put 7 Z = 1 + ir. Condition (1) ensures that we may choose a positive number 
e z , such that claim (a) in (ii) is satisfied. Indeed, otherwise we could choose a sequence 
(t n ) of positive numbers, such that t n — > as n — > oo, and v a (x) + t n r = \m(z + t n j z ) < 
v a (Re(z + t n j z )) = v a (x + t n ) for all n. This implies that 

v a (x) + t n r - v a (x) v a (x + t n ) - v a (x) 
r = < , 

for all n, which contradicts (1) and the fact that the right hand side converges to v' a (x) 
as n — > oo. 

Regarding assertion (b) in (ii), we remark first that statement (ii) in Lemma 13.11 
ensures that H a (z + tj z ) G C + for all t in (0, e z ). We note next that 

lm(^ H a (z + t lz )) =\ m (HUzH) = ((\ m (H' a (z)),Re(H' a (z))),(l,r)). 



I t=0 



Condition (2) thus ensures that Im^i^^ + tj z )) ^ at t = 0, which implies (b). This 
completes the proof. ■ 

3.6 Lemma. Let a be a strictly positive number, let z be a point in S a \ {— c a }, and put 
£ = H a (z) = z + azG^(z) E R (cf. Lemma\J^iii)). 

Then the Cauchy transform G Ua of v a has the non-tangential limit - at More 
precisely for any positive number S we have that 

lim G Ua {£ + w) = \ 

w— >0 y 

where As is given by (12. 7p . 

Proof. By Lemma [3To1 we may choose 7 2 in C and e z in (0, oo), such that z + t^ z e 
for all t in (0, e z ), and such that H a (z + tj z ) —> H a (z) = £ non-tangentially (from C + to 
R) as 1 1 0. Using Proposition I3.2lf ii) it follows that 

lim G Va (H a (z + tj z )) = lim ■ 



40 Q ' " " 40 z + 17 2 z 



(note that z / 0, since ^ 9a)- Applying then Lemma 12.11 (to the function w 
—G Va (£ + w)), we may conclude that actually 

lim G Va (( + w) = - 

w— f g 

for any positive number 5, as desired. ■ 

For any a in (0, oo) we introduce next the function P a : R — > IR (cf. Lemma l3.lT iii)) 
given by 

P a (x) = H a (x + iv a (x)), (x e R). (3.14) 

In particular we put 

s a = P a {-c a ) = H a (-c a ). (3.15) 

In the following lemma we collect some properties of P a that will be needed in the sequel. 
We defer the rather technical proof to Appendix [A] 
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3.7 Lemma. For any positive number a the function P a : R — > R has the following 
properties: 

(i) P a is continuous on R and analytic onl \ {— c a }. 

(ii) P a satisfies that 



x 



+ « + «/o° S dt ' 



if x < —c n 



2x + a-a J °° ( ,_,gi;^^ 2 dt, if x > -c a . 



(3.16) 



(x-t) 2 +i; ci (x) 

(iii) The number s a := P a {— c a ) is strictly positive. 

(iv) lim (x + a — P a (x)) = 0. 

X— >oo 

(v) P a is a strictly increasing bisection of R onto R, and P' a (x) > for all x in R\ {— c a }. 




The graphs of the functions P\/2, Pi, P2 and P 
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3.8 Theorem. Let a be a positive number, and consider the function Q a : 
defined by 

Qa(x) = , Va{X) (X G ' 



£ + v a {x) 



R ->• [0,oo) 
(3.17) 



We then iave 



(i) The free Gamma distribution v a with parameter a is absolutely continuous (with 
respect to Lebesgue measure). 

(ii) The density f a of u a is given by 



/«(0 



0, if f G (-00, s a ] 

iQ^Pi- 1 ^)), if £e(s«, 00), 



(3.18) 



where Pa ^ denotes the inverse of P a (cf. (v) in Lemma \3.7\) . 
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(iii) The support of v a is [s a , oo). 

(iv) f a is analytic on (s a , oo). 




The graphs of the densities of the free gamma distributions V\/2 , v\ and, respectively, 1/2, via- 

Proof of Theorem lff.#L (i) As described in Subsection 12.3} the singular part of v a 
(with respect to Lebesgue measure) is concentrated on the set 

S = {£eR|lim^ |O^ + h/)| = oo}. 

From Lemma [3.71 it follows that P a is a continuous, increasing bijection of R onto R. By 
(13.141) this implies that H a maps Sa bijectively onto R. For any £ in R \ {s a } it follows 
then by Lemma 13.61 that 

MmG Va ^ + iy) = -, (3.19) 

2/4-0 Z 

where z is the unique point on Sa \ {— c a }, such that H a (z) = £. We may therefore 
conclude that S C {s a }, and the proof of (i) is completed, if we verify that u a has 
no atom at s a . As mentioned in Subsection 12.11 it follows from (the proof of) [BV93t 
Proposition 5.12] that v a has at most one atom, which, if it exists, is necessarily equal to 
the non-tangential limit of the Voiculescu transform <fi Ua {z) as z — > 0, z G C + . Note here 
for z in C + that by (12.51) and Proposition I3.2( i). 

<Pv a ( z ) = z ^v a (5) = azG^ (z) = a + a / dt 

Jo z — t 

(cf. formula (13. ip ). Hence, by dominated convergence, 

/ r°° te _< \ 
lim u) Vn (iy) = lim [a + a dt ) = a — a = 0. 

It follows that the only possible atom for u a is 0, and since s a > (cf. Lemma [3.7( h)). 
we conclude that v a has no singular part. 
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(ii) By Stieltjes inversion (cf. Subsection I2.3p . the formula 

/ a (0 = --limlm(GW(£ + h/)), 

7T 2/4-0 

produces an almost everywhere defined density for v a with respect to Lebesgue measure. 
According to (13.191) we have for all £ in M. \ {s a } that 

/«(e) = --im(-), 

7T \zJ 

where z is the unique point in Sa \ {— c a }, such that H a (z) = £. Writing z = x + iv a (x) 
for some (unique) x in R\ {— c a }, we have that 

£ = H a (x + iv a (x)) = P a (x), so that x = Pj^^iOi 

and therefore 

fa(0 = --lm( 1 ) = -^^L = ig a(x) = iQ^Pi- 1 )^)). 

The proof of (ii) is completed by noting that if £ < s a , then x < — c a , so that f a (x) = 0, 
and therefore /«(£) = by the previous calculation. 

(iii) This is an immediate consequence of (i), (ii) and the fact that v a (and hence Q a ) 
is strictly positive on (— c a ,oo) (cf. Lemma ETTT iv)). 

(iv) Since v a is analytic on (— c a , oo) (cf. Lemma ETffl i)). it follows immediately from 
(I3.17P that so is Q a . By (i) and (v) of Lemma I3T7I the function P a is analytic on (—c a , oo) 
with strictly positive derivative. This implies that Pa ^ is analytic on (s a ,oo), and 
altogether we thus conclude that f a = -Q a o P^ l) is analytic on (s a ,oo). ■ 



4 Behavior at the limits of the support 

In this section we study the behavior of the density f a of v a around the lower bound s a 
of the support and at infinity. We start with the latter aspect. 

4.1 Proposition. Let a be a positive number, and consider the density f a of v a (cf. 
fl3~THD ). We then have 

v /a(0 a 

hm — — - = ae . 

Proof. Consider the function P a introduced in (13.141) . Since P a is a strictly increasing 
bijection of R onto M., it suffices to prove that 

lim f a (P a (x))P a (x)e p ^ =ae a . 



Using Theorem I3.8f ii). Lemma l3.4l (iv) and Lemma l3.7f iv). we find that 

1 Vg(x) 
x^oo 71 X 2 + V a (x) 



hm f a (P a (x))P a (x)e p ^ = lim - ^ P Q (x)e p ^> 



lim V ^ \f" a W]f^W| c -^W 



— ■ 1 ■ air • 1 ■ e = ae , 

7T 



as desired. 
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We turn next to the behavior of / Q (£) as £ 4- s a . We study initially how s a varies as a 
function of a. 

4.2 Proposition. For any positive number a consider the function H a and the quantities 
c a and s a defined by (l3.ip . fl3.4B and (13.151) . respectively We then have 

(i) H a satisfies the differential equation: 

H' a {z) = a + z + (z' 1 - l)H a (z), (zeC\ [0, oo)). 

(iii) H>>(-C a ) = 1 - f < 0. 

(iv) lim Q ^ s a = 0, and lim^oo s a = oo. 

(v) c a is an analytic function of a, and ^ = ^-^-c 2 ) • 

(vi) s a is an analytic function of a, and ^ = . In particular s a is a strictly 
increasing function of a . 

Proof, (i) Differentiation in the first equality in (13. ip and partial integration leads to 

oo e -t ,00 e -t 



H'Jz) = l + a I -dt-az — dt 

z — t J [z — t) 



( r e~* i 00 f 

z~ 1 H a (z) — az[ + / 

Uz-tJo J 



-oo -t 

dt 



' ^ : - t 

z~ l H a (z) + a - (H a (z) - z) = a + z + (z^ 1 - l)H a (z) 



for all z in C \ [0, oo). 

(ii) From Lemma 13. 3\ (i) and ( I3.15P it follows that 

= H' a {-c a ) = a-c a + (-c" 1 - l)H a (-c a ) = a - c a - ±±^s c 

from which (i) follows immediately. 

(iii) Differentiation in (i) leads to the formula: 

H';(z) = 1 - z- 2 H a (z) + (z- 1 - l)H' a (z), (zeC\ [0, oo)). 
Combining this with Lemma 13. 3^ we find that 

K(-c a ) = 1 - c- 2 H a (-c a ) + = 1- c~ 2 s a . 
At the same time it follows from ( 13. 2 p that 



= *»/„ r _^d. = -2«y (^*<0, 



and thus (iii) is established. 

(iv) From Lemma IBTlT i) we know that c a — > as a — > 0, and that c a — > oo as a — > oo. 
In combination with (ii) and (iii), respectively, it follows that s a has the same properties. 



15 



(v) For any x in (0, oo) it follows from (I3.2p -( 13T3]) that 



H' a {-x) = l-af J^-- 2 dt = l-aF(-x) 



so that 



F{-x) = ~(l-H' a {-x)), and F'{-x) = --H'^(-x), (ie(0,oo)). (4.1) 
a a 

Using (iii) and (ii) we find thus that 

F'(-Ca) = --K{-Ca) = -(0- 2 S a - 1) = -{ " " °\ - 1^ - " " ' ' 1 



a a aVcQ,(l + c a ) / aca,(l + c a ) 

In particular we see from (iii) that F'(—c a ) > 0, and (14.11) shows that F is analytic on 
(—oo,0). From the defining formula: F(— c a ) = - and the implicit function theorem (for 
analytic functions; see |FG02t Theorem 7.6]) it follows thus that c a is an analytic function 
of a with derivative given by 

dcQ, 1 c Q (l -\- Cq,) ^_ 00)) 

da a 2 F'(—c a ) a(a — 2c a — c 2 ) ' 

(vi) From (ii) and (v) it is clear that s a is an analytic function of a. Consider now 
the function S: (0, 00) x (0, 00) — > K given by 

c 

5(c, a) = (a — c), (c, a G (0, 00)), 

and note that s a = S(c a , a) for all positive a, and that 
dS a — 2c — c 2 c 

gc" (c ' a)= (i + c) 2 ' and ^ (c ' a) = TT^' ( c ' ae (°'°°))- 

Using the chain rule and (v) it follows thus that 

ds a dS , .dc a dS , s a — 2c a — c 2 n c a (l + c a ) c a 
(c a ,a)- h7-(Ca,a) = — - ; - — + 



da dc da da (1 + c a ) 2 a(a — 2c a — c 2 ) 1 + c D 
c a c a c a (a "I - 1) 



a(l + c a ) 1 + c a a(l + c a ) 
and this completes the proof. ■ 

Let a be a real number contained in an interval /. For functions g, h: I — > C, such 
that ^ h(I), we use in the following proposition the notation: u g(x) ~ h(x) as 2 — > a" 
to express that lim x _;. a |M = 1 . 

4.3 Proposition. For any positive number a we put 

6#«(-Ca) 



70 ^'(-^ 



Tii en 7a > 0, and we have that 
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(i) v a (x) ~ 7a 12 \x + Cq) 1 / 2 , as x J, — c a . 

(ii) lim P " (X ]~ g " = -\laK{-c a ) > 0. 

(iii) / Q (0 ^=^(£ - s a ) 1/2 , as £ 4 s a . 

Proof. Using formula (13. 2p we note first for any A; in {2, 3, 4,...} that 

r°° f e -t 

Hi k \~c a ) = -ak\ / - +1 dt < 0, (4.2) 

and in particular this verifies that 7 a > 0. 

(i) Using formula (IA.4[) (in Appendix |2J we find that 

^K(x) 2 ) = 2v a (x)v^(x) — > ° — = la, 

dx xi ~ Ca Jo W^ dt 

where the last equality results from (14. 2p . Since v a (—c a ) = 0, it follows from the above 
calculation and the mean value theorem that v Q (x) 2 ~ 7 a (x + c a ) as x j. — c a , and this 
proves (i). 

(ii) Using that H a (—c a ) = s a and H a (—c a ) = (cf. Lemma l3~3|) we find by Taylor 
expansion that 

P a (x) = Re(H a (x + iv a (x))) = Re(s a + ^H'^(—c a )(x+c a +iv a (x)) 2 +o{jx+c a +iv a (x)\ 2 ) 
and hence by application of (i), 

p a (x) ~ s a 1 / v a (x) 2 \ o((x + c a f + v a {xf ) . 

= ~H a (-c a )(x + c a — J + > -2# a (-c Q )7 a , 

as x | — c a . Since 7 Q > 0, formula (14. 2p shows that the resulting expression above is 
positive, and hence (ii) is established. 

(iii) Recall from Theorem 13.81 that 

f a (P a (x)) = - X (X E [~C a , OO)). 

7T x i + v a {xy 

By application of (i) it follows thus that 

1/2 



f a {P a {x)) ~ -^-(x + c a ) 1/2 , as 1 1 -c a . (4.3) 
Using (ii) we have also that 



and hence 



2 

P i~ 1} (0 + Ca~ 777/7 r(f-Sa), &S f | -S a . 



Combining this with (14. 3 p we find that 

f (ft 7a /2 (i J i" 1> (0+C tt ) 1/2 ,1/2 



< 7r^(-^(-0) 1/a 



as £ 4- s a- Applying finally Proposition I4.2f ni). we obtain (iii). 
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5 Unimodality 



In this section we establish unimodality of the densities f a . We start with a few prepara- 
tory results. 

5.1 Lemma. For each positive number R, let (0, it) — > (0, oo) be the function given 
by 

$ R (0) = F(Rsm(0) e w ), (0 e (0,tt)), 

where F is the function introduced in (I3.3P . Then for any R in (0, oo) there exists a unique 
number Or in (0, tc) such that is strictly decreasing on (0, 6r] and strictly increasing 
on [0 R , it). 

Proof. We note first that for any r in (0, oo) and in (— tt, it] we have, using the change 
of variables t = ru, that 

F(re w )= !°° — dt= [°° rwr™ 

V ; J (rcos(0) -t) 2 + r 2 sm 2 (0) J r 2 (cos(0) - u) 2 + r 2 sin 2 (0) 

du. 



/ 1 - 2u cos(0) + u 2 
Hence, for a fixed positive number R we have that 

f°° ,, lp -uRsm(6) 

1 — 2u cos(#) + u 
Then define the function ^/ R : (—1, 1) — >■ (0, oo) by 



®r(0) = F(Rsm(0)e w ) = \ - , ~ 2 du, (0 e (0,vr)). 



roo -uRVls 2 

so that $r(0) = ^r(cos(0)) for in (0, tt). Since the function h-> cos(#) is strictly 
decreasing on (0, tc), it suffices then to show that ^r is strictly decreasing on (—1, t)r] and 
strictly increasing on [tjr, 1) for some number rjR in ( — 1, 1). For this we consider for any 
u in (0, oo) the function ipR ;U : (—1, 1) — > (0, oo) given by 



By a standard application of the theorem on differentiation under the integral sign, it 
follows that ^r is differentiable on (—1,1) with derivative 

r°° d 

^r( s )=J o ^R,u(s)du, (se(-l,l)). (5.1) 
For any u in (0, oo) and s in (—1, 1) we note further that 



£ HMs)) = £ ( ln(u) - uRVT^ - ln(l - 2us + u 2 )) 
= uRs(l - s 2 y 1/2 + 2u(l - 2ms + m 2 )- 1 , 
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so that 



ds 2 



ln(^fl„(s)) = uR(l - s 2 )- 1 ' 2 + uRs 2 {l - s 2 )~ 3/2 + Au 2 {l - 2us + u 2 \ 



= uR(l - s 2 )" 3/2 + 4u 2 (l - 2us + u 2 Y 2 > 0. 

Since 



\n(ip Ru {s)) - 



ds 2 ' ^R,u( S ) 1>R,M 2 

we may thus conclude that ip Ru > and hence that ip' Ru is strictly increasing on (—1,1). 
Since this holds for all u in (0, oo), it follows further from (15. ip that ^>' R is strictly increasing 
on ( — 1, 1). Thus, ty R is either strictly increasing, strictly decreasing or of the form asserted 
above. However, by Fatou's Lemma, 

f°° u f°° u 
liminf ^ R (s) > / —du = oo, and liminf \&r(s) > / — du = oo, 

sti J Q (1 - u) 2 4-1 J (1 + u) 2 

and hence ty R must have the claimed form. ■ 

5.2 Lemma. Let a be a strictly positive number, and consider the functions Q a , P a and 
f a given in flXTTj) . fl3TT5j) and (ETTgjl . We then have 

(i) For any p in (0, oo) the equation: 

Q a {x) = p 

has at most two solutions in (— c a , oo). 

(ii) For any p in (0, oo) the equation: 

= p 

has at most two solutions in (s a , oo). 

Proof, (i) Let p be a strictly positive number, and assume that there exist three distinct 
points xi,x 2 ,x 3 in (— c a ,oo) such that 

p = Qa(xj) = -lm( 1 ), (j = 1,2,3). 

It is elementary to check that the points z in C \ {0}, for which — lm(l/z) = p, consti- 
tute the circle C p in C with center ^-i and radius ^- (except for the origin). Thus our 
assumption is that C p intersects the set 9' Q (given in ( I3.1UP ) at three distinct points. Note 
that 

C p = U(i + ^) I P e *}} = (i(cos(/3) + i(l + sin(/3)) | /3 G (-vr, vr]}. 

Writing a point ^(cos(/3) + i(l + sin(/3)) from C p \ {0} in polar coordinates re l8 (r > 0, 
9 G (0, 7r)), it follows that 

rsin(#) = i(l+sin(/3)), and r 2 = ^( C os 2 (/3) + l+sin 2 (/3)+2 sin(/3)) = ^(l+ s in(/3)), 
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so that 

l + sin(/3) rsin(fl) 1 . 
r = = = - sin S . 

2/rr pr p 

Hence, 

C p ={±sin(#)e ie |flG(0,7r]}, 

and our assumption thus implies that there are three distinct points 9i, 62,63 in (0, it), 
such that - sin(6'j)e l6,J G 9' a , j = 1,2,3. According to (I3.10p and (13. 7p . this means that 
the equation 

F(±sin(#)e ie ) = - (5.2) 

has (at least) three distinct solutions in (0,7r). However, Lemma 15.11 asserts that the 
function 

$ Vp (0) = F(±sin(0)e i0 ), (0G(O,7r)), 

is strictly decreasing on (0, and strictly increasing on [9y p , ir) for some 9y p in (0, 7r). 
Hence the equation (15.21) has at most two solutions in (0,7r), and we have reached the 
desired contradiction. 

(ii) Let p be a strictly positive number, and assume that there exist three distinct 
£1, £2, £3 in (s a , 00) such that f a (£j) — P> j — 1; 2, 3. Then there exist three distinct points 
^1,^2,^3 in (— c a , 00), such that P a (a;j) = £j> 3 — 1>2, 3, and it follows from formula 
(13381) that 

P = fa{Pa{x )) = ^QaiXj), (j = 1, 2, 3). 

This contradicts (i), and the proof is completed. ■ 

5.3 Theorem. For each a in (0, oo) the density f a of the free Gamma distribution v a is 
unimodal. In fact, there exists a number co a in (s a , oo) such that f a is strictly increasing 
on [s a ,uj a ] and strictly decreasing on [u Q , oo). 

Proof. The proof is an elementary consequence of Lemma f5.2f ii). but for completeness 
we provide the details. We know that that f a is continuous, that /<*(£) > whenever 
£ > s a , and that 

f a (s a ) = = lim / a (£) 

(cf. Lemma l3.1l (iv). Theorem 13.81 and Proposition 14. ip . In particular it follows that f a 
attains a strictly positive maximum at some point u a in (s a , oo). We show next that f a 
is non- decreasing on [s a , cu a ). Indeed, if this was not the case, we could choose £1,^2 in 
(s a ,u) a ) such that 

6 < and > / Q (£ 2 ). 

Choosing an arbitrary number p in (f a (^2), /a(£i)), it follows then by continuity of f a 
that there must exist si in (s a , £1), S2 in (£1,^2) and S3 in (^2,^0) such that 

fa(si) = p, (i = 1,2,3). 

Since this contradicts Lemma l5\2T ii). we conclude that f a is non-decreasing on [s^u^]. 
This further implies that f a is strictly increasing on that same interval, since otherwise f a 
would be constant on a non-degenerate sub- interval, which is precluded by Lemma loT2T ii). 

Similar (symmetric) arguments show that f a is strictly decreasing on [u a , 00), and this 
completes the proof. ■ 
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6 Asymptotic behavior as a — > 

In this section we study the asymptotic behavior of the free Gamma distributions u a , as 
a 1 0. We start by considering convergence in moments. 

6.1 Proposition. The measures converge in moments to the measure t _1 e _< l(o i0 o) (t) dt 
as a 10. More precisely we have for any p in N that 



1 

a 



t p u a (dt) — > / t p_1 e -t dt, as a 10. 



Proof. It follows from Proposition 14. ll that v a has moments of all orders (cf. also [BG06J). 
It follows moreover from [AnO!) Lemma 6.5] that for all p in N the free cumulant r p (a) 
of v a equals the classical cumulant c p (a) of fi a (the classical Gamma distribution with 
parameter a). The latter cumulants may be identified by considering the Taylor expansion 
at of \og(fi a (u)) . Using dominated convergence, it follows that for any u in (—1, 1) we 
have that (cf. ( 1532)) ) 



log(£«(«)) = a ^ (e iui - 1) °— dt = a f°° ( V ±!^!_) e -* dt = a V 

from which we may deduce that 

r p (a) = Cp(ot) = a(p — 1)! for all p in N. 

Using the Moment-Cumulant Formula (cf. |NiSp| ) it follows further that the p'th moment 
rrip(a) of v a is given by 



k=2 ^ ' <n ,...,<Jt>l 



r qi {a)r q2 {a)---r qk {a) 



k=2 v 7 g 1 ,...,9 fe >i 

91 H H9fc=P 



for any p in N. In particular we see that m p (a) is a polynomial in a of degree p with no 
constant term and linear term a(p — 1)!. For any p in N we may thus conclude that 

/oo /»oo 
V v a {dt) = lm p {a) — >(p-l)! = j( t^Vdt, 

as desired. ■ 

We show next that the densities of ^z/ a actually converge point- wise to £ -1 e - *l(o )0 o)(£) 
as a j. 0. 

6.2 Lemma. Consider the functions P a defined in (13.141) . 

(i) For any x in (0, oo) we have that P a (x) — > x, as a 1 0. 

(ii) For any 7/ in (0, oo) we nave that Pi _1 '(y) ?/, as a | 0. 

Proof, (i) Let x be a fixed number in (0, oo). From (13. ip and (I3.14p it follows that 

r°° te~ l 

P a (x) = x + w a (x) + a + a / — dt, (a e (0,oo)). 

x — t + w a (x) 
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Lemma l3.4f iii) clearly implies that v a (x) — > 0, as a — > 0, and hence it suffices to show 
that 

dt — >0, as a -> 0. (6.1) 



x — t + w a {x) 

From Lemma l3.4( iii) it follows furthermore that we may choose ot\ in (0, oo), such that 
> f xe~ x , whenever a G (0, at]]. Then for all t in (0, oo) and a in (0, a%] we have that 



a 



< — — — < 



x — t + iv a (x) v a (x) I a 
For any t in (0, oo) \ {x} we note further that 



(6.2) 



a 



x — £ + if a (x) 



< a-, 



\x — 1\ 



0, as a —7- 0. 



(6.3) 



Combining ( 16. 2 p and ( 16. 3 p it follows by dominated convergence that ( 16. ip holds, as desired. 

(ii) Let y in (0, oo) and e in (0, y) be given. From (i) we know that P a (y — e) y — e, 
and P a (y + e) — )■ y + e, as a — )• 0. Hence we may choose «2 in (0, oo) such that 

-Pa(y — e) < 2/, an d -P«(y + e) > whenever a G (0, a?\. 

Then for any a in (0, ctz] we have that 

ye [P a (y-e),P a (y + e)} = P a ([y - e,y + e}), (6.4) 

since P a is increasing and continuous. It follows from (16 ,4p that 

P^ 1 \y) G [y — e,y + e], whenever a G (0, a 2 ], 

and since e was chosen arbitrarily in (0,y), this establishes (ii). ■ 

6.3 Proposition. For any x in (0, oo) we have that 

-f a (%) ~> x~ l e~ x ', astt->0. 

Proof. Let x be a fixed number in (0,oo), and note that Lemma l3.4( iii) implies that 
Va ^ — > irxe~ x , as a — » 0. Using then formula (I3.18P we find that 



J/a(^a(x)) 



v a (x)/a 



xe 



tt(x 2 + f a (x) 2 ) x 2 + 
It suffices thus to show that 



as a — > 0. 



- / Q (P a (x)) - / a (x) 
a 1 1 



0, as a — > 0. 



For all positive a we put y a := Pi (a?) 5 and Lemma l6T27 ii) then asserts that y a — >■ x, as 
a — >• 0. Given any number 5 in (0, x) we may thus choose a\ in (0, oo) such that 



y a G [5, oo), whenever a G (0, cci]. 



(6.5) 
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For any a in (0, a±] we find then by application of (I3.18P that 



-\fa(Pa(x)) -f a (x)\ = -\f a (P a (x)) - f a (P a (y a ))\ 

1 cr 1 



a 



< 



v a (x)/a v a (y a )/a 



x 2 + v a (x) 2 yl + v a (y a y 
v a (x)/a - v a (y a )/a | v a {y 



(6.6) 



x 2 + v a (x) 2 



x 2 + v a (x) 2 yl + v a (y a y 



Consider now in addition an arbitrary number 7 in (0,1). By Lemma l3.4( iii) we may 
then choose a<i in (0,«i], such that 



sup 

m£[i5,oo) 



v a (u) 



— nue 



a 



< 7, whenever a G (0, a 2 ] 



(6.7) 



Using (16. 5 p and (16. 7p we find that 

— — — < Ky a e~ Va + 7 < 7r sup ue~ u + 1 < 00, whenever a G (0, a 2 ]- 

a «e(o,oo) 

Together with the fact that y a — > x as a — > 0, this implies that 

1 _ 1 1 agcw0 

y 2 a+V a (y a ) 2 ~ y2 a + a 2^^lY 2 ^' ^ ^ 

Since also x2+ J — ► js, as a — > 0, another application of ( 16. 8 p then yields that 



(6.8) 



v a {y 



a 



x 2 + v a (x) 2 y 2 a + v a (y a ) 2 



0, as a — > 0. 



In view of (16. 6p it remains thus to show that 

v a (x)/a - v a (y a )/a 



0, as a —7- 0. 



(6.9) 



x 2 + v a (x) 2 

For this, note that for any a in (0, 02] we have by new applications of ( 16 .5p and ( 16. 7p that 



v a (x)/a - v a (y a )/a 



x 2 + v a (xy 



< 



27 + 7r|xe x — y a e Va \ 



x' 



Since u t- >■ ue u is continuous at x, we may choose in (0, a^] such that \xe x — y a e Va \ < 
7r _1 7, whenever a G (0,a 3 ], and then 



v a (x)/a - v a (y a )/a 



x 2 + v a (x) 2 



37 

< — , whenever a G (0, a 3 ]. 



Since 7 was chosen arbitrarily in (0, 1), this verifies (16. 9p and completes the proof. 
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A Proofs of Lemmas 13.41 and 13.7 



In this appendix we provide detailed (but rather technical) proofs of Lemma 13.41 and 
Lemma [3.71 We start with the following preparatory result: 

A.l Lemma. Let a be a positive number and consider the function v a : R — >■ [0, oo) 
given by fl3.6p - fl3.7p . We then have 



(i) If < e < x, then 



Vry \ X . 



> 2a (x — e)e x e arctan( — 



(ii) For any e in (0, 1) we have for all sufficiently large x that 

V a (x) < 2a ( X + € ) e -x+e arctan (_ 



1 — e \v a (x)J' 
Proof, (i) Recall first (cf. (J32D) that 



a J (x - t) 2 + v a (x) 



2 dt, (x G (-c a , oo)), (A.l) 



Assume next that < e < x. Using flA.ip and the change of variables u = we find 

that 

1 f x+t te"' , (x-e)e- x - € f x+e 1 
> / -, ^ ^ dt > ^ /— — / , . . dt 



a (t-a;) 2 + ^(x) 2 ^(x) 2 1 + (^fy) 2 

(A.2) 

t/Va(x) 1 , N , 2(x - e)e~ x - e , , , 

t> Q (x) du = — — arctan( 



WaW 2 J-e/v a (x) 1 + ^ Ua(^) Ua(a:) ' 

from which the desired estimate follows immediately. 

(ii) Let e be a given number in (0,1), and note that for any t in (0, oo) and x in 

(e, oo), 

te"* , 4 te~* 



l[o,»-e](0j 77 vf7 r^l[ic+e,oo)(£) < e te . 



(t-x) 2 +^(x) 2 l^"^^' (t_ x )2 +VQ ( x )i 

Hence, by dominated convergence, 



to"* , f°° te-* 

dt, / — ; — ttt dt — > 0, as X —7- oo. 



/ (t - x) 2 + v a (x) 2 J x+e (t - x) 2 + v a (x) 
Thus, for all sufficiently large x we have by 0A.10 that 

1_ , te~* (z + e)e- x+£ f x+e 1 

( _e) «-y,- £ (t-x) 2 +^(x) 2 * - Wq (x) 2 x_ e i+(^) 2 * 

Q( j (A.3) 

= 2(X t(x) arCtan( ^) } ' 
which yields the desired estimate. ■ 
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Proof of Lemma \3-4[ (i) Consider the function F: 1 x (0, oo) — > M. given by 

f°° te~* 

F{x,y) = F{x + iy) = dt, (x G M, y > 0). 

Vo (a - t) 2 + r 

Using formula (13. Sp in the case a = 1, it follows that 

= 1 -^Im^x + k/)), Glx (0,oo)), 

and since the imaginary part of an analytic function is analytic (as a function of two real 
variables), we may conclude from this that F is analytic onlx (0, oo). By differentiation 
under the integral sign we find in particular that 

aF r°° (x - t)te- f 

— {x,y)--2J o {{x _ t)2 + y2)2 dt, 



and 



9F, , n f°° te~* 

— {x,y) = -2y dt. 

<9y Jo ((a; - t) 2 + y 2 ) 2 



Since u a (:r) > and F(x,v a (x)) = ^ for all x in (— c a ,oo), and since ^(x,y) < for 
all (x,y) in R x (0, oo), it follows then from the implicit function theorem (for analytic 
functions; see [FG02I Theorem 7.6]), that v a is analytic on (—c a , oo) with derivative given 
by 



dF ( I \\ _ f°° (x-t)te 1 i, 

<{*) = ~f {x ' Va{x)) = : j vr <)2+ ^ )2)2 „ (x e ( -c Q , oo)). (a.4) 



%{x,v a {x)) v a (x)f™ df 



In particular v a is continuous on (— c a ,oo). From the defining relations ( lA.lj) and f!3.4p 
it is standard to check that v a (x) — > as x 4- — c a . Thus, t> Q is continuous at — c a as well 
and hence on all of M. 

(ii) Using Lemma IA.l( i) we find for any positive e that 

v ( r ) 2a(x — e)e~ x ~ e arctan(— f-?) 
hmmf > hmmf L - L - = 2ae — = e air, 

x-^oo XG x x^foo XG x 2 

where we have used that v a (x) — > as x — > oo (cf. Lemma [A.l( ii)). Letting then e — > 0, 
it follows that 

liminf > avr. (A.5) 

x— >-oo xe 31 

Using Lemma [A .1 H i) we find similarly for e in (0, 1) that 

7 > (V) 2a(x + e]G~ x+t arctan(— %-^) P e 

hmsup — — < hmsup — — L - L - = air, 

x— >oo xg x x^>oo (1 e)xe x 1 e 

and letting then e — > 0, we conclude that 

limsup fy^l < a;7r. (A. 6) 

Combining (IA.5[) and (IA.6[) completes the proof of (ii). 
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(iii) Let e be a fixed number in (0, |]. For any x in [e, oo), we note then that 



te 



o (x - t) 2 + v a (x)' 



dt + 



tc 



x+e 



(x — t) 2 + v a (x)'- 



dt 



) 



px—e poo poo 

< ae~ 2 ( / tc~ l dt + / £e~* dt) < ae~ 2 / te~ t dt = ae 

^ JO J x+e ' JO 



and thus 



xd \e,oo 



iP ) a Wo (x - t) 2 + v a (x) 2 



dt 



tc- 



x+e 



(x — t) 2 + v a (x)' 



dt 



0, as a —)• 0. 



In combination with (lA.lj) this implies that we may choose ot\ in (0, oo) such that for all 
a in (0, «i] and all x in [e, oo) we have that 



1 - e < a 



tc 



dt 



, (x - t) 2 + v a (x) 2 
= 2a(x + e)e~ x+€ v a (xy 1 arctan (— frr) < Tra(x + e)e~ x+t v a (xy 1 , 

^ Hot \p£ ) ' 

where we have re-used the calculation (1A.3j) . Hence, it follows that 
v a (x) tt(x + e)e~ x+e 



< 



for all x in [e, oo) and a in (0, cti]. 



a 1 — e 

Since e < |, we find in particular for all a in (0, a±] that 

sup v a (x) < K e a, where K e := 2iiy/e sup (x + \)C X < oo. 

xd[e,oo) xS[e,oo) 

For any x in [e, oo) and a in (0, oo), we note next that 

px+e te - t 



(A.7) 



(A.8) 



1 > a 



dt > a(x — e)e x t v n (x) 



x+e 



dt 



, (x - i) 2 + v a (x) 2 

= 2a(x - e)e~ x ~ e v a (xy 1 arctan (^y). 

In combination with fl A.8j) this shows that for all x in [e, oo) and a in (0, cti] we have that 

v a (x) 



a 



> 2(x - e)e x e arctan(eif e a 



i„.-i^ 



Hence, we may choose «2 in (0, ai], such that 
v a (x) 



a 



> 2(x - e)e- x - € {\ - e)f = (1 - e)ii(x - e)e~ x - 



(A.9) 



for all x in [e, oo) and a in (0, ol^[. Combining now flA.7j) and (IA.9|) . it follows for any a 
in (0, a<2\ that 



sup 

x£[e,oo) 



v a (x) 



ITXC 



a 



< 7T SUp 

x£[e,oo) 



(xe~ x - (1 - e){x - e)e~ x - t ) V ((1 - e)~ x {x + e) 



e x+e -xc~ x ) 



(A.10) 
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Using that the function x h- > xe x is bounded on (0, oo), it is standard to check that 
sup (xe~ x -{l-e)(x~e)e- x - t ), sup ((1 - e)" 1 ^ + e)e~ x+e - xe~ x ) — ► 0. (All) 

xe(0,oo) x6(0,oo) 



To complete the proof, assume that positive numbers 5 and 7 are given. By (lA.llj) we 
may then choose e in (0, 5 A |] such that the right hand side of flA.lOp is smaller than 7. 
Applying the above considerations to this e, it follows that we may choose a 2 in (0, 00), 
such that 

^ S Lip 

xe[e,oo) 



sup 

x£[5,oo) 



— 7rxe 



a 



nxe x 



a 



<7, 



whenever a G (0,0:2]. 

For the proof of Lemma [3.71 we need the following auxiliary result. 
A. 2 Lemma. Let a, r, e be positive numbers such that e < 1. We then have 

(i) x -rj- ^dt — >0, as x ^00. 

(ii) If < e < x, then 



2(x-e) r e- x - e , , f x+e Te"* 



f Q (x) (ar -t) 2 + v a (x) 

2(x + e) r e- x+e 
< — arctan 



(iii) For all sufficiently large positive x we have that 



ax 2 (x 2 + e) Jx-i- (x — t) 2 + v a (x) 2 ax[x 2 — e 

Proof, (i) We note first that for a; in, say, (2, 00), we have that 



2x t r e- f , , _ „ « /" 2x 1 



dt < (2x) r e~ x -^ I ji^dt = 2 r e~ 2 x r (x 3 - ex)e 



L (x-t) 2 + v a (x) 2 " v ' J x+ * m2 

X X 

and similarly that 



dt < (x - t-) r e- x ' 2 / dt = e" 2 (x - f ) r (f - ex)e-^ 2 . 



, /2 ( x -t)2 + Va ( x )2 " v x > J x/2 ^ 
Moreover, 

dt < / — — dt = — / t r e~ r dt, 



2x ( X ^) 2 ~t~ ^a(^) 2 J2x a ' 2 ^ 



2,r 



and 

,x/2 ^ r ,x/2 t r e -r 4 /-x/2 

dt < / — — dt = — / t r e~ r dt 



(X - t) 2 + V a (x) 2 -Jo (f)2 X2 7 Q 

Now, the sum of the left hand sides of the 4 estimates above is equal to the integral 
f, n , e i ei 7 — tt^ — dt, and the sum of the right hand sides is clearly of size oix^ 1 ) 
as x — > oo. This shows (i). 
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(ii) Assume that < e < x. Arguing as in the proof of Lemma l3.4l (ii). we find that 

re"* , (x + e) r e- x+t f x+e 1 

dt < - Hr- / , , , dt 



(t - x y + Va ( x y ~" - Va {xf J x _ e 1 + (^fj) 2 

2(x + e) r e- x+ < e 
= — arctanl — HrJi 

which proves the second estimate in (ii). The first estimate follows similarly. 

(iii) Considering x in (l,oo), we find by application of (ii) and Lemma lA.ll fi) (with e 
replaced by ejx) that 



t 2 e"* . 2(x + ^) 2 e~ x ^ 



dt < x —- arctan( : 



(x — t) 2 + V a (x) 2 V a (x) xv a (x) 

2(x+f) 2 e-^arctan(^ R ) + c )* e ' 



< 



2a(x — -)e x i arctan( — Vr) ax(x 2 — e) 



which proves the second estimate in (iii). Regarding the first estimate, we note first that 
it follows from (i) that 

tor* , e 

dt< — 



! {0 ,oo)\(x-± ,x+± ) (x - t) 2 + v a (x) 2 ax 
for all sufficiently large x, and hence by (lA.lj) and (ii) 

< / 7 — ttt dt < ^— r arctanl 



a (x — t) 2 + v a (x) 2 v a (x) xv a (x)< 

X 

for all sufficiently large x. For such x we may thus conclude that 

, . 2a(x + ^)e- x+ -* 
v a {x) < : — e arctan( 



\ — l \xv a (x) ' 



which in combination with (ii) yields that 



t 2 e"* , 2(x - ^) 2 e- x -^ 

■ dt > x —— arctanl 



(x — t) 2 + V a (x) 2 V a (x) ""^xv a (x) 

2(1 -~)(x- ^) 2 e~ x -^ arctan 



~~ 2a(x + -)e _a;+ ^ arctan( 



- xv a (x) ' 

(x — e)(x 2 — e) 2 e - ^ 



ax 2 (x 2 + e) 

for all sufficiently large x. This completes the proof. ■ 

Proof of Lemma \3. 7[ (i) Since H a is analytic on C \ [0, oo) and v a is analytic on 
R \ { — c a }, it follows immediately from (I3.14p that P a is analytic onl\ {— c a }. Since v a 
is continuous on M, it follows also that so is P a . 
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(ii) For x in (—00, —c a ), formula (13.1 6p follows immediately from (13. ip . since v a (x) = 
0. For x in [—c a , 00) we find, using Lemma l3~TT (iii). (13. ip and ( lA.ip . that 

P a (x) = H a (x + iv a (x)) = Re(H a (x + iv a (x))) 

x + a + a I Re ( ^— — \ dt 

-x + w a (x) — t/ 

f°° (x - t)te~* 
= x + a + a / — --^ dt 

j {x - ty + v a {x) 2 

f 00 te~ l , r t 2 e~* 
= x + a + ax / ? -— - — - dt-a -— - dt 

j {x-ty + v a (x) 2 Jo {x-ty + v a {xy 

= 2x + a-a / — — dt, 

Jo [x - t) 1 + v a {x) 2 

as desired. 

(iii) Considering the function H a restricted to (— 00, 0), it follows from (13. ip and 
dominated convergence that 

lim H a (z) = —00, and lim H a (z) = 0. 

2— 00 2—^0 

26K 2S(-oo,0) 

From (13. 2p it follows further that 

H' a {z) > 0, if z G (-00, -Ca), and H' a (z) < 0, if 2 G (-c„, 0). 

Since P Q = ifa, on (— 00,— c a ], we deduce from these observations that P a is strictly 
increasing on (—00, — c a ], and that s a > inf ze (_ Caj0 ) H a (z) = 0. 

(iv) Using formula (I3.16P as well as (i) and (iii) of Lemma IA.2I we find for any e in 
(0, 1) that 

/ f°° t 2 e~* 
lim sup (x + a — P a (x)) = lim sup ( — x + a / -— ^ dt 

= hmsup^-x + a / (x _ t)2 + v (x] 2 dt 

x 

„2 1 J\2, " 



< lim sup I — x + a — — 

x->do v ax(x 2 - e) J 

= limsup(-x + ^(l + ^)e^) 

= hmsup Ue* -l)+(^- + l+ ^fcy)et) 



- 1 

e « — 1 

lim — : = 2e. 
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Arguing similarly we find next that 

(x - e)(x 2 - e) 2 e~ 



liminf (x + a — P a (x)) > liminf ( — x + 

x— >-oo x— >oo \ 



a- 



ax 2 (x 2 + e) 
= l^(-* + <^(l- =& ) e -* 

= liminf f x(e"^ - 1) - -^T* + e2 ~ x r x * £ ( 1 - -^rle"- 
= -3e. 

Combining the two estimates obtained above, and letting e — > 0, we obtain (i). 

(v) From (13.1 6p and (iv), it is clear that P a (x) — > ±oo as x — > ±oo, and since P a is 
continuous, it suffices thus to prove that P' a (x) > for all x in R \ {— c a }. In the proof of 
(iii) we already noted that P' a (x) > for all x in (— oo, —c a ). For x in (—c a , oo) we find 
by differentiation in (13.14j) that 

P^x) = H' a (x + iv a (x))(l + iv , a (x)) 

(A.12) 

= Re(H' a (x + iv a (x))) -\m(H' a (x + iv a (x)))v' a (x), 

where we have used that P' a {x) G R, so that 

= lm(i^(z)) = Re(^(x + iv a (x)))v' a (x) + \m(H' a (x + iv a (x))). (A.13) 

According to Lemma 13731 Re(H' a (x + iv a (x))) > 0, and hence (1A.13I) implies that 

-\m(H' a (x + iv a (x))) 
Va[X) Re{H' a (x + iv a (x))) ' 

which inserted into flA.12[) yields that 

P> (x) = Re(H> (x + iv (x))) + HH' a (x + iv«(x))) 2 = \H> a (x + iv a (x))\> 
a{ ) [ a{ + M ))) + Re(H' a (x + iv a (x))) Re(H> a (x + iv a (x))) > °' 

as desired. ■ 
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